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Abstract. 

We consider the propagation of galactic cosmic rays under 
assumption that the interstellar medium is a fractal one. An 
anomalous diffusion equation in terms of fractional deriva- 
tives is used to describe of cosmic ray propagation. The 
anomaly in used model results from large free paths ("Levy 
flights") of particles between galactic inhomogeneities and 
long time a particle stays in a trap. An asymptotical solution 
of the anomalous diffusion equation for point instantaneous 
and impulse sources with inverse power spectrum relating to 
supernova bursts is found. It covers both subdiffusive and su- 
perdiffusive regimes and is expressed in terms of the stable 
distributions. The energy dependence of spectral exponent of 
observed particles in different regimes is discussed. 



1 Introduction 

A key to understanding the mechanism for cosmic-ray ori- 
gin and acceleration is determination how cosmic-ray parti- 
cles propagate through the interstellar medium (ISM). If the 
propagation process is determined by scattering at magnetic 
field inhomogeneities which have small-scale characters and 
can be considered as homogeneous Poisson ensemble, it can 
be described by a normal diffusion model (Ginzburg and Sy- 
rovatskii (1964); Berezinsky et al. (1990)). The normal dif- 
fusion process is characterised by a mean-squared displace- 
ment that increases with time, (r^ (t)} ex t, and by a Gaussian 
propagator. 

However, during a few last decades, many evidences from 
both theory and observations of the existence of multiscale 
structures in the Galaxy have been found (see, for exam- 
ple, Lee and Jokipii (1976); Kaplan and Pikelner (1979); 
Lozinskaya (1986); Ruzmaikin et al. (1988); Vanshtein et 
al. (1989); Bochkarev (1990); Falgarone et al. (1991); 
Burlaga et al. (1993); Molchanov et al. (1993); Meish and 
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Bally (1994); Armstrong et al. (1995); Minter and Spangler 
(1996); Cadavid et al. (1999)). Filaments, shefls, clouds are 
entities widely spread in the ISM. A rich variety of structures 
that are created in interacting phases having different proper- 
ties can be related to the fundamental property of turbulence 
called intermittency. 

The stretching, bending and folding of magnetic field lines 
by turbulent motions of a medium partially coupled to the 
field make magnetic field also highly intermittent, especially 
on the smaller scales. As the turbulent zones do not fill 
space, the tool of the fractal geometry to characterise the ISM 
should be used. In such a fractal-like ISM we certainly do not 
expect the normal diffusion to hold. 

Generahsation of this equation leads to anomalous diffu- 
sion (see the reviews by Bouchaud and Georges (1990); 
Isichenko (1992); West and Deering (1994) and Uchaikin 
and Zolotarev (1999)). In the case of anomalous diffusion 
(r^ (t)) (xt^,t^ oo, where the exponent 7 differs from 1, a 
value, that corresponds to the normal diffusion. Anomalous 
diffusion has a non-Gaussian propagator. 

In our recent papers (Lagutin et al. (2000, 2001a,b)) we 
proposed an anomalous diffusion (superdiffusion) model for 
solution of the "knee" problem in primary cosmic-rays spec- 
trum and explanation of different values of spectral exponent 
of protons and other nuclei at £' ^^ 10^ -^ 10^ GeV/nucleon. 
The anomaly results from large free paths ("Levy flights") of 
particles between magnetic domains — traps of the returned 
type. These paths are distributed according to inverse power 
law oc Ar~^~°', r — + 00, a < 2 being an intrinsic property 
of fractal structures. We also assumed that the mean time a 
particle stays in a trap, (r), was finite. 

In this paper we consider the propagation of galactic cos- 
mic rays in fractal ISM without the latter assumption. We 
suppose now that the particle can spend anomalously a long 
time in the trap. An anomalously long time means that (r) = 
L dr T qij) = 00, so the distribution of the particles stay- 
ing in traps, (/(r), has a tail of power-law type oc Bt~^~^, 
i — > 00 with /3 < 1 ("Levy trapping time"). 



2 Particle diffusion in a fractal medium 

After acceleration in a source the particle can be in one of 
two states: a state of "Levy flights" or a state of rest-state of 
motion in a trap. Diffusion in this model is a process in which 
the particle state changes successively at random moments in 
time. Based on the continuous time random walk theory of 
MontroU and Weiss (1965) or on integral equation (see, for 
example, Uchaikin and Zolotarev (1999); Uchaikin (1999)) 
we can readily derive the fractional diffusion equation. The 
equation for Green's function G(r, t, E; Eq) without energy 
losses and nuclear interactions under condition that the par- 
ticle started from origin Tq = at time to = with energy 
Eq has the form 



dG 

dt 



-D{E, a, l3)Dl-\-/\r/^G{r, t, E; E^) 



-S{r)S{t)S{E - Eo) 



(1) 



Here Dq^ denotes the Riemann-Liouville fractional deriva- 



tive (Samko et al. (1987)) 



DUit) 



1 d 
r(l - ij.)di 



(t-T)-^f{T)dT, ^X<l, 



(—A)"' — the fractional Laplacian (called "Riss" operator 
(Samko etal. (1987))) 



(-A)"/^/(x) = -i 



Ai/(x) 



dy, 



where l> a,x ^ R™, y e R™, 



fe=0 



and 



drra{^)= j{l-e^y)'\v\- 

R™ 



'dy. 



The anomalous diffusivity D{E, oi, P) is determined by the 
constants A and B in the asymptotic behaviour for "Levy 
flights" (A) and "Levy waiting time" (B) distributions: 

D{E,a,l3) (xA{a)/B{E,P). 

Taking into account that the probability to stay in a trap 
during the time interval t for particle with charge Z and mass 
number A depends on particle rigidity as oc R^^, we find 
D oc {E/Zy. 

The solutions of equations (1) with zero boundary condi- 
tions at infinity can be found by Laplace-Fourier transforma- 
tions with use of formulae (Samko et al. (1987)) 



^-\tp,ij. 



D^^Gir,t,E;Eo)dt 



A^ / e-^'G{r, t, E; Eq) dt = A^G(r, A, E; Eq), 
Jo 



I e*'='-(-A)"/2G(r, t, E; Eq) dr 



R3 



= |fc|" / e'^''G{r, t, E; Eq) dr = \k\°'G{k, t, E; Eq). 

R3 

The Laplace-Fourier transformation solution of (1) is 

G(fc, A, E] Eq) = 6{E - Eo)\'^~^ 

OO 

X I exp{-[Xf^ + D{E,a,(3)\k\"]y)dy. 



Invert Laplace-Fourier transform we find 

G{r,E,t;EQ) ^ S{E ^ EQ){D{EQ,a,l3)t''y^'' 

x^^t'^(\r\{DiEo,a,P)t^)-'^''), (2) 

where 

OO 

*("'« (r) - f q*") (rr^) q[^-^^ (r) t^^^'^ dr. (3) 



Here q"^' (r) — (27r)^^ J exp(— i/cr — \k\")dk is the den- 
sity of three-dimensional spherically-symmetrical stable dis- 
tribution with characteristic exponent a < 2 (Zolotarev et 
al. (1999); Uchaikin and Zolotarev (1999)) and qf'^^ (t) is 
one-sided stable distribution with characteristic exponent (3 
(Zolotarev (1983)): 



Qi 



f '^^ (t) = (27rz)-i / exp(Ai - X^)dX 



.(2) 



Let us remember that qj^ (r) is the normal (Gaussian) dis- 
tribution density, q^ ' (r) is the three-dimensional Cauchy 
density [n{l + r^)]^^, qi (r) is Levy-Smirnov distri- 

bution. Other stable densities cannot be expressed through 
elementary functions, but there exist representations in terms 
of convergent and asymptotic series (Uchaikin and Zolotarev 
(1999)). Based on equation for Green's function (1) it is easy 
to formulate the fractional diffusion equation for concentra- 
tion: 

^ = -D{E,a,l3)Dl-^{-A)"^^N{r,t,E) 



+S{r,t,E), 



(4) 



where S(r, t, E) is a density function of sources distribution. 



3 Spectra 



4 Conclusion 



Using Green's function (2) we can find the cosmic ray con- 
centration for sources interesting for astrophysics. Thus, for 
point instantaneous source with inverse power spectrum 

S{r,t,E) = SQE-PS(r)S{t) 

we have 

N{r,t,E) ^ SoE-P{D{E,a,P)t'^y^^" 

x^i'^'^^{r{DiE,a,P)tP)-'^"). (5) 

For point impulse source 

S{r, t, E) = S'oE-PS(r)e(T - t)e{t) 

^(^) ^ 1 0, r < 0, 
the solution is of the form 



Nir,t,E) 



S'qE^p 



D{E,a,f3)r/» 



-3/3/0 



max[0, t-T] 

x¥,-'^U\r\{DiE,a,(3yy'/''' 



(6) 



Using the representation TV — NqE ^ and the asymptotic 
behaviour of the scaling function vI/("''^) (r) (a < 2, f3 < 1) 

,(",/3) 



*3 (^^0)ocr-(3-")^ 



^3 ' (^ ^ °°) °^ ''" 



-3-a 



one can evaluate the variation of spectral exponent Arj — 

ri{E ^ cx)) - r]{E -^ 0). It follows from (5) that 

N{r,t,E)(xE-P+^, E^O, 

N{r,t,E) (xE-P-\ E-^oo, 

A77 = 2(5. 

In other words, in case of point instantaneous source in both 
subdiffusive (/3 < a/2) andsuperdiffusive(/3 > q;/2) regimes 
the spectral exponent of observed particle increases with en- 
ergy on 2(5, i.e. the cosmic ray spectrum steepens (the "knee"). 
The similar estimates for the impulse source (see (6)) give 

N{r,t,E)(xE-P+^, E^O, 

N{r,t,E)oiE-P-^'^, E^oo, 

A77 = (5(l + l//3). 

Our analytical and numerical studies however show that this 
property of energy spectrum (the "knee") is lacking in the 
regimes of normal diffusion (a = 2, (3 = 1) and subdiffusion 

(a = 2, /3 < 1). 



We considered the propagation of galactic cosmic ray in the 
fractal interstellar medium. Anomalous diffusion equation in 
terms of fractional derivates describing of cosmic ray pro- 
pagation has been formulated. An asymptotical solution of 
this equation covered both subdiffusive and superdiffusive 
regimes has been expressed in terms of stable distributions. 

Our results showed that the "knee" in the primary cosmic 
ray spectrum is due to anomalously large free paths ("Levy 
flights") of particles, being an intrinsic property of the fractal 
interstellar medium. 
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